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Abstract. The STAR experiment at the Relativistic Heavy-Ion ColUder at Brookhaven 
National Laboratory is carrying out a spin physics program in high-energy polarized p -\- p 
collisions at a/s = 200 — 500 GeV to gain a deeper insight into the spin structure and dynamics 
of the proton. 

One of the main objectives of the spin physics program at RHIC is the extraction of the 
polarized gluon distribution function based on measurements of gluon initiated processes, such 
as hadron and jet production. The STAR detector is well suited for the reconstruction of various 
final states involving jets, 7r°, tt^, e^ and 7, which allows to measure several different processes. 
Recent results will be shown on the measurement of jet production and hadron production at 
^ = 200 GeV. 

The RHIC spin physics program has recently completed the first data taking period in 2009 
of polarized p+p collisions at = 500 GeV. This opens a new era in the study of the spin- flavor 
structure of the proton based on the production of VK"*^^-* bosons. Recent STAR results on the 
first measurement of W boson production in polarized p-\-p collisions will be shown. 



1. Introduction 

The spin structure and dynamics of the nucleon is one of the fundamental and unresolved 
questions in Quantum Chromodynamics (QCD). Contrary to atomic physics where the total 
angular momentum of an atom can be accounted for by basic quantum theory in terms of its 
underlying atomic constituents, a complete theoretical description of the proton spin decomposed 
in terms of contributions from quark and gluon angular momenta and spins is still missing. 
Various experimental programs have been conducted in the past to deepen our understanding 
on the proton spin. Deep-inelastic scattering (DIS) experiments have clearly established that 
the quark spin contribution is small and accounts for only ^ 25% of the proton spin [HIS]. 

High energy polarized p + p collisions at ^/s = 200 — 500 GeV at RHIC provide a unique 
way to probe the proton spin structure and dynamics using hard scattering processes [3]. 
The production of jets and hadrons is the prime focus of the gluon polarization studies. The 
production of VF~(+) bosons at ^/s = 500 GeV provides an ideal tool to study the spin-flavor 
structure of the proton. This has been pointed out at the very early design stages of the polarized 
proton collider program at RHIC 14] . 

The first global analysis of polarized DIS data, as well as results obtained by the PHENIX 
[5] and STAR [6J experiments in polarized p^p collisions at RHIC placed a strong constraint 
on the gluon spin contribution in the gluon momentum range of 0.05 < x < 0.2, and suggested 
that the gluon spin contribution is not large in that range [7J . Constraining the polarized gluon 
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Figure 1. Recent STAR 2006 All inclusive jet result. 



distribution function, Ag, through inclusive measurements has been, so far, the prime focus of 
the STAR physics analysis program of the Run 3/4 [8j, Run 5 [6J and Run 6 [9j data samples. 
Inclusive measurements, such as inclusive jet production, integrate over a fairly large x region for 
a given jet transverse momentum region. While those measurements provide a strong constraint 
on the value of Ag integrated over a range in x, those measurements do not permit a direct 
sensitivity to the actual x dependence. This motivates the need for correlation measurements 
in polarized p + p collisions. 

The STAR collaboration has presented a first measurement of the longitudinal spin transfer 
Dll in inclusive A (A ^ P^~) and A (A ^ pn^) production in polarized proton-proton collisions 
at a center-of-mass energy of ^/s = 200 GeV [lOj. This measurement may provide constraints 
on strange (anti) quark polarization pjj and can yield new insight into polarized fragmentation 
functions. 

The first data taking period in 2009 of polarized p^p collisions at ^/s = 500 GeV opens 
a new era in the study of the spin-flavor structure of the proton based on the production of 
bosons. W~^^'^ bosons are produced predominantly through u + d {u + d) collisions 
and can be detected through their leptonic decay. Quark and anti-quark helicity distribution 
functions are probed at large scales set by the mass of the W boson (Q ^ mw) where theoretical 
calculations are well under control. The theoretical framework is well developed to describe the 
production of W bosons in high-energy polarized p + p collisions including the description of 
the W decay leptons subject to experimental cuts, such as the transverse momentum pT and 
pseudo-rapidity 77 of the final-state leptons [121 [13] . This development profits from a rich history 
of hadroproduction of weak bosons at the CERN SPS and the FNAL Tevatron and provides a 
firm basis to use W production as a new high-energy probe in polarized p^p collisions [14j. 

The STAR Collaboration has provided critical measurements on the sudy of transverse spin 
effects in polarized p+ p collisions at RHIC. An overview of these results is discussed in a 
separate paper [15]. 
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Figure 2. Recent STAR 2006 All neutral pion result (a) and STAR 2006 All results for 
charged pion (ti~ and tt^) / jet correlation measurements (h/c). 



2. Recent jet All results 

STAR reconstructs jets with the midpoint cone algorithm using clusters of charged track 
momenta measured with the STAR Time Projection Chamber (TPC) and tower energy 
deposits in the STAR Barrel Electromagnetic Calorimeter (BEMC) within a cone radius 
of i? = \/ /S.T]'^ + A(/)2 [9J. The jet sample was required to be within a fiducial range of 
—0.7 < //Detector < 0.9 for the 2006 data sample with a cone radius of i? = 0.7. The dominant 
fraction of jet events in the 2006 data sample are based on a jet patch ( JP) trigger that required 
a minimum energy deposition for a group of towers over a region of A77 x A(/) = 1.0 x 1.0. This 
trigger was taken in coincidence with a minimum-bias condition using the STAR Beam-Beam 
Counter (BBC). 

Throughout the following discussion, four gluon polarization scenarios have been used as 
input to NLO perturbative QCD calculations of All- The GRSV standard case refers to 
a global analysis fit of polarized DIS data [16]. The case for a vanishing gluon polarization 
(GRSV-ZERO) and the case of a maximally positive (GRSV-MAX) or negative (GRSV-MIN) 
gluon polarization have also been considered. The GS-C [17j set of is reflected by a large 
positive gluon polarization at low x, a node around x ^ 0.1 and a negative gluon polarization 
at large x at a scale oi ^ 1 GeV^. 

Figure [1] shows the most recent STAR preliminary result of All for inclusive jet production 
as a function of pT based on the 2006 data sample of 4.7 pb~^. Typical beam polarization values 
during the 2006 data-taking period were 55 — 60%. The leading systematic uncertainty comes 
from the trigger and jet reconstruction biases. Differences between observed and true jet pT are 
estimated using PYTHIA and GEANT simulations and result in corrections being applied to the 
jet pt values. Additionally, the trigger conditions at STAR can artificially bias the analyzed data 
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Figure 3. Comparison of STAR 2005 di-jet data and a PYTHIA Monte- Carlo showing the 
yield at the top and the ratio at the bottom for three characteristic di-jet variables (See text for 
further details). 



sample towards particular flavors of partonic collisions. A conservative systematic uncertainty is 
evaluated to account for this effect which incorporates all allowable gluon polarization scenarios. 
The All curves in Figure [T] are derived from NLO fits to world polarized DIS data by two 
separate theory groups, GRSV [16j and GS [17]. The currently measured result is consistent 
with a small gluon polarization scenario for the STAR x-range of 0.03 < x < 0.3. This x-range 
accounts for - 50% of the total AG integral, i.e. AG = f^ Agdx, for GRSV-STD. 

3. Recent neutral and charged pion All results 

Neutral pions are identified via the two-photon decay channel, which accounts for more than 
98% of TT^ decays. Two primary sub-detectors, the STAR Barrel Electromagnetic Calorimeter 
(BEMC) and respective Shower Maximum Detector (SMD) are employed to reconstruct decay 
photons from mid rapidity neutral pions. Charged pion production is of particular interest 
since the difference of the longitudinal double-spin asymmetries for tt^ and tt" production, 
- tracks the sign of A^, due to the opposite signs of the polarized distribution 

functions for up and down quarks. The STAR Time Projection Chamber (TPC) offers robust 
reconstruction and identification of charged pions up to a transverse momentum of 15 GeV/c. 
Particle identification is accomplished using measurements of ionization energy loss in the TPC. 

Figure [2] (a) shows All for neutral pion production at mid rapidity as a function of pT 
based on the 2006 data sample |18j. This new measurement is shown in comparison to NLO 
calculations assuming different gluon polarization scenarios [9l [18] . This measurement reaches 
to higher pt values compared to previous measurements by PHENIX [19] and STAR [20j. At the 
present level of precision, the data excludes extreme gluon polarization scenarios. The STAR 
collaboration has recently presented also results for neutral pion production at forward rapidity 
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Figure 4. (a) Statistical precision of the longitudinal double-spin asymmetry, All, cls a function 
of the di-jet invariant mass, M, for the 2006 RHIC data sample, (h) Longitudinal double-spin 
asymmetry. All, for di-jet production as function of the ratio M/y^ for different topological 
combinations of the STAR BEMC and the STAR EEMC acceptance region. 

based on the STAR Electromagnetic Endcap Calorimeter (EEMC) and the STAR Forward Pion 
Detector (FPD) [21j. These measurements provide an important milestone for future photon-jet 
coincidence measurements at STAR. 

Figure [2] (b) and (c) show STAR's new preliminary result of All for charged pion (tt" and 
7r+) /jet correlation measurements at mid rapidity. Charged pions are reconstructed opposite 
a jet that triggered the experiment. The data are shown versus z = Pt{^) /prij^t) and were 
obtained in 2006 [22]. Full NLO calculations for this observable have recently been released 
[23] . These calculations will provide the basis for charged pion results to be included in a global 
analysis. The data shown for charged pion production are compared to a LO MC evaluation of 
All excluding extreme scenarios of A^. The measurement of is of particular interest, 

since its analyzing power is large because of the large ug scattering contribution to the production 
cross section and sizable Au/u. In the future, more precise measurements of All(^^) have a 
great potential for providing a better understanding of Ag. 

4. Status and prospects of di-jet production 

Correlation measurements such as those for di-jet production allow for a better constraint of 
the partonic kinematics and thus the shape of Ag. At LO, the di-jet invariant mass, M, is 
proportional to the product of the x values of the partons, M = v^y^^i^' whereas the pseudo- 
rapidity sum of the final-state jets, 773 + 774, is proportional to the logarithm of the ratio of the 
X values, 773 + 774 = ln(xi/x2). Photon-jet coincidence measurements are expected to provide 
a theoretically clean way to extract A^. A LO extraction of alone would allow a model- 
independent way to constrain the x dependence, which would be an important contribution, 
without making an a priori assumption on the functional form of A^ as is currently required 
in a global analysis. Measurements at both ^/s = 200 GeV and ^/s = 500 GeV are preferred to 
maximize the kinematic reach in x and possibly provide a means to observe effects of scaling 
violations at fixed x by measuring different pT values. The wide acceptance of the STAR 
experiment permits reconstruction of di-jet events with different topological configurations. 
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Figure 5. Event displays from the 2009 data taking period at ^/s = 500 GeV showing a W 
candidate event (a), a candidate event and a di-jet candidate event (c). 

i.e. different 773/774 combinations, ranging from symmetric (xi = X2) partonic collisions to 
asymmetric (xi < X2 or xi > X2) partonic collisions. This, together with the variation of the 
center-of-mass energy, is expected to constrain Ag over a wide range in x of approximately 
^ 2 • 10~^ < X < 0.3 for di-jet and photon-jet events. The NLO framework for correlation 
measurements exists and therefore those measurements can be used in a global analysis [24j. 

Figure 3 shows a comparison of STAR 2005 di-jet data and a PYTHIA Monte-Carlo 
simulation with the yield at the top and the ratio at the bottom for the invariant mass M, the 
pseudo-rapidity sum 773 + 774 and the cosine of the partonic center-of-mass scattering angle cos ^* 
determined from the pseudo-rapidities of both final-state jets. The normalization between data 
and Monte-Carlo is fixed by the invariant mass distributions. The same normalization factor is 
then applied to both other di-jet distributions. Data and Monte-Carlo are found to be in good 
agreement . 

Figure 4 (a) shows the statistical precision of the longitudinal double-spin asymmetry. All, 
as a function of the di-jet invariant mass, M. Those uncertainties, extracted from the 2006 data 
sample, are compared to a LO MC evaluation of All computed with a PYTHIA MC sample 
using different event weights to account for different polarized gluon distribution functions of 
GRSV ^ and DSSV 0- The size of the statistical uncertainty at the highest invariant mass 
bin is at the level of the difference between GRSV-STD and DSSV. 

Figure 4 (b) shows the expected precision for the longitudinal double-spin asymmetry, All^ 
for di-jet production as a function of Mj ^ for different topological combinations of the STAR 
BEMC and the STAR Endcap Electromagnetic Calorimeter (EEMC) acceptance regions. Taking 
into account the different 77 ranges covered, and equivalently, the different cos0* regions being 
probed, each panel represents a different range in Xi/x2. At LO, cos ^* amounts to tanh ^ ^3-yy4 ^ 
The upper left panel effectively probes asymmetric partonic collisions where predominantly a 
low-x gluon collides with a high-x quark at large invariant masses. The effective variation of 
xi and X2 amounts to 0.2 < xi < 0.6 and 0.07 < X2 < 0.2. In contrast, a kinematic region of 
larger x values is probed at predominantly symmetric partonic collisions such as the one shown 
in the lower right panel. The effective variation of xi and X2 is roughly equal and given by the 
horizontal axis of the lower right panel. The projected uncertainties are shown for a luminosity 
of 50 pb~^ and a beam polarization of 60%. Those projected uncertainties are compared to a LO 
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Figure 6. Lego plot showing the BEMC Et and the TPC pt distribution for a W candidate 
event (a) and a di-jet candidate event (h). 

evaluation of All and a full NLO All calculation. Scale uncertainties are shown as a shaded 
band for DSSV and GRSV-STD reflecting a variation of the invariant mass M as a hard scale 
of 2M and 0.5M. Asymmetric cuts are imposed for the LO MC and the NLO determination of 
min(pT) > 7GeV/c and max(pT) > lOGeV/c. The result of a LO MC evaluation using GS-C 
[TT] for is also shown. GS-C is still consistent with the current inclusive jet results [9J. A 
cone radius of i? = 0.7 has been used. Good agreement is found between a LO MC evaluation 
of All and a full NLO calculation. Scale uncertainties are found to be small in comparison to 
the variation of the chosen polarized gluon distribution functions, in particular, at large values 
of M. 

While inclusive measurements from STAR have provided important constraints on A^, and 
will continue to do so, the impact of correlation measurements to the x dependence should 
greatly enhance our understanding of A^. 

5. Measurement of W boson production at STAR 

The data used for the first W boson production analysis at STAR were collected in 2009 colliding 
polarized proton beams at 250 GeV. The STAR detector systems used in this analysis are 
the STAR Time Projection Chamber (TPC), the STAR Barrel Electromagnetic Calorimeter 




Figure 7. Illustration of the W offline selection criteria (top) and the distribution of the TPC 
Pt (a) the BEMC Et (h) and both BEMC SMD layers (c) / (d) in azimuth and pseudo-rapidity. 



(BEMC) and STAR Electromagnetic Endcap Calorimeter (EEMC). Only their tower response 
has been taken into account in this analysis. The BEMC was used to measure the transverse 
energy, Et, of e^. The suppression of the QCD background by several orders of magnitude was 
based on the TPC, BEMC and EEMC. Figure 5 shows event displays from the 2009 data taking 
period for a W candidate event (a), a candidate event and a di-jet candidate event (c). A 
lego plot is shown in Figure 6 displaying the BEMC Et and the TPC pt distribution for a W 
candidate event at the top and a di-jet candidate event at the bottom. 

Proton-proton collision events focusing on W production at ^/s — 500 were identified by a 
two-step energy requirement in the BEMC. Electrons and positrons from W production at mid 
rapidity are characterized by large Et peaked at ^ Mw/'^. (Jacobian peak). At the hardware 
trigger level (LO), a high tower (HT) calorimetric trigger condition required Et > 7.3 GeV 
in a single BEMC tower. At the software trigger level (L2), a dedicated trigger algorithm 
was developed that required that a 2 x 2 tower cluster Et sum exceeds 13 GeV. The offline 
selection of W candidate events is based on kinematical and topological difference between 
leptonic decays, and QCD background events as shown in Figure 5 and Figure 6. Events 
from + u decays contain a nearly isolated with a neutrino in the opposite direction 

in azimuth. The neutrino escapes detection, leading to a large missing energy. 

Figure 7 illustrates the basic idea of the selection of W candidate events. The 2D histograms 
show for one W candidate event the TPC pt (a) and the BEMC tower Et (b) distribution 
along with both BEMC shower- maximum detector (SMD) layers (c) / (d) as a function of the 
azimuthal angle and the pseudo-rapidity 77. The approximate area of a 2 x 2 tower cluster is 
indicated by a small circle around the large energy deposition in a single BEMC tower. Also 
shown is the size of a cone (near-cone) around the electron candidate. The away-side area 
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Figure 8. Evolution ofW offline selection cuts (See text for further details). 



is sketched by the black dashed rectangular region. The TPC high-p^^ entry appears in the 
same r] — (/) region as the large BEMC Et entry. Both BEMC SMD layers show a large energy 
deposition at the same location in the rj — (j) plane as the TPC high-p^^ and the BEMC high- 
Et entry. An electron candidate is defined to be any TPC track with pT > lOGeV/c that is 
associated with a primary vertex with \z\ < 100cm, where z is the distance along the beam 
direction. A 2 x 2 BEMC tower cluster Et sum is required to be larger than 15 GeV. The excess 
BEMC Et sum in a 4 x 4 tower cluster centered around the respective 2x2 tower cluster is 
required to be below 5%. In addition, the distance between the 2x2 cluster tower centroid 
and the TPC track is required to be less than 7 cm. A near-cone is formed around the electron 
candidate direction with a radius in Tj-cj) space of i? = 0.7. The excess BEMC and TPC Et sum 
is required to be less than 12% of the 2x2 cluster Et- The away side requirement is based 
on a cut on the BEMC and TPC Et sum to be less than 8 GeV. This sum is extended over 
the full 77 range of the BEMC and EEMC and Ac/) = 0.7 as shown in Figure 7. Figure 8 shows 
the evolution of all cuts for the BEMC 2x2 Et distribution. A clear Jacobian peak emerges 
characteristic for W production in contrast to QCD background dominating the low Et region. 

A data-driven procedure was applied to estimate the level of QCD background events. 
The first part consists of an estimate of the impact of the missing calorimeter coverage for 
— 2<77<— Ion the cluster Et distribution. This estimate was obtained by performing a 
parallel analysis without the EEMC as an active detector. The difference in yield of accepted 
events as a function of Et provides an absolute estimate for the missing calorimetric coverage. 
The remaining background is estimated by fitting a separate QCD shape distribution in Et to the 
region of Et < 19 GeV shown as the black histogram in Figure 9. The fitted QCD background 
distribution is shown by the blue histogram in Figure 9. The background subtracted distribution 
is shown by the yellow histogram. Figure 9 also shows the distribution from a PYTHIA Monte- 
Carlo simulation. Good agreement in yield and shape is found between the data and Monte- 
Carlo distribution for the charge-sum BEMC Et distribution. In the meantime, a considerable 
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Figure 9. Reconstructed BEMC lepton transverse energy for W events. 



effort has been performed on the TPC cahbration to allow a full charge-sign separation. This 
and in particular the first measurement of the respective charge-separated cross-sections and 
parity- violating singe-spin asymmetries Al s^t mid rapidity will be discussed in a forthcoming 
publication. 

Future planned high-statistics measurements at forward and backward pseudo-rapidities 
will focus on constraining the polarization of d and u quarks, respectively to deepen our 
understanding of the spin structure of the QCD sea. 
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